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We use silicon strip detectors (originally developed for the CLEO III high energy particle physics 
experiment) to measure fluid particle trajectories in turbulence with temporal resolution of up to 
70,000 frames per second. This high frame rate allows the Kolmogorov time scale of a turbulent 
water flow to be fully resolved for 140 > R\ > 970. Particle trajectories exhibiting accelerations 
up to 16,000 ms^'^ (40 times the rms value) are routinely observed. The probability density 
function of the acceleration is found to have Reynolds number dependent stretched exponential 
tails. The moments of the acceleration distribution are calculated. The scaling of the acceleration 
component variance with the energy dissipation is found to be consistent with the results for 
low Reynolds number direct numerical simulations, and with the K41 based Heisenberg-Yaglom 
prediction for R\ > 500. The acceleration flatness is found to increase with Reynolds number, 
and to exceed 60 at Bx = 970. The coupling of the acceleration to the large scale anisotropy is 
found to be large at low Reynolds number and to decrease as the Reynolds number increases, 
but to persist at all Reynolds numbers measured. The dependence of the acceleration variance 
on the size and density of the tracer particles is measured. The autocorrelation function of an 
acceleration component is measured, and is found to scale with the Kolmogorov time t^. 



1. Introduction 

Fluid turbulence may be characterized in terms of variables defined at points fixed in space 
(the Eulerian reference frame), or in terms of the trajectories of fluid particles (the Lagrangian 
reference frame). This distinction applies both to theoretical formulations of turbulence and to 
experimental techniques for characterizing turbulent flows. Although the formulation of fluid 
dynamics is generally considered to be more tractable in terms of Eulerian variables, the critical 
issues of transport and mixing in turbulence are more directly related to the properties of fluid 
trajectories (Shraiman & Siggia 2000; Yeung 2001/?; Sawford 2001), and are often addressed us- 
ing Lagrangian techniques (Pope 1994). There are also many applications in which the transport 
or aggregation of particulate matter in turbulence is important in its own right, such as water 
droplet aggregation in clouds (Vaillancourt & Yau 2000) or the industrial production of nanopar- 
ticles (Pratsinis & Srinivas 1996). 

In a basic sense, data obtained from Lagrangian and Eulerian measurements are complemen- 
tary. In the Eulerian frame, one is typically concerned with differences between quantities (ve- 
locity component, scalar concentration, etc.) measured for several points separated by a fixed 
distance in space. It would be equally interesting to study time differences, however such mea- 
surements made within the Eulerian framework are difficult to interpret, since the large scales 
of the flow will sweep the turbulence past a fixed detector, causing temporal structure to be en- 
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tangled with spatial structure. The true temporal structure of turbulence is only revealed when 
fluctuations are measured along a particle trajectory, in the Lagrangian frame. 

Turbulence has traditionally been studied in the Eulerian frame rather than the Lagrangian 
frame for technical reasons. The hot wire anemometer used in conjunction with the Taylor 
frozen flow hypothesis provides extremely accurate Eulerian data in turbulent gas flows over 
a broad range of Reynolds number, but no comparably effective technique has been available 
for Lagrangian measurements. However, as the range of Reynolds number accessible to direct 
numerical simulations (DNS) of turbulence has expanded, numerical studies of particle trajec- 
tories in isotropic turbulence have yielded important insights (Vedula & Yeung 1999; Yeung 
2001a; Gotoh & Rogallo 1999; Gotoh & Fukayama 2001). Although very little experimental 
data for Lagrangian properties in fully developed turbulence has been available, this is beginning 
to change. Basic issues, such as the Richardson law for particle dispersion (Moniti & Yaglom 
1975), the Heisenberg-Yaglom prediction of fluid particle accelerations (Heisenberg 1948; Ya- 
glom 1949), and Kolmogorov scaling of temporal velocity differences have remained untested 
for many decades. Recent experiments by Ott & Mann (2000), Voth et al. (1998), and Mordant 
et al. (2001) have made significant progress in addressing each of these but there are still ma- 
jor hmitations. In particular, hmitations in spatial and temporal measurement resolution have 
obscured the small scales at large Reynolds number. More complete information about particle 
trajectories in real (non-idealized) turbulent flows is needed to guide the development of models 
of transport in applications (Pope 1994). 

In principle, fluid particle trajectories are easily measured by seeding a turbulent flow with 
minute tracer particles and following their three-dimensional motions with an imaging system. 
In practice, this is a very challenging task because changes in particle velocity or acceleration 
can take place on time scales of order of the Kolmogorov time, — (i^/e)^/^ where v is the 
kinematic viscosity and e is the energy dissipation per unit mass. In order to observe universal 
scaling behaviour we require that the Reynolds number, defined by Re = uL/v, where u is the 
rms velocity and L is the energy injection scale, approach 10^. (This is equivalent to requiring 
that the Taylor microscale Reynolds number R\ = (15Re)^/^ approach 1000.) In a laboratory 
water flow {v « 10~^ m^ s~^) with convenient energy injection scale (L « 0. 1 m) and assuming 
e = u^/L we must have e w 10 m^ s~^, which impUes r w 0.3 ms. In order to measure particle 
accelerations, motions which take place over this time scale must be fully resolved. 

Although conventional imaging systems based on charge coupled devices (CCD) have been 
used for three-dimensional tracking of particles in low Reynolds number flows (Snyder & Lum- 
ley 1971; Sato & Yamamoto 1987; Virant & Dracos 1997; Ott & Mann 2000) they do not provide 
adequate temporal resolution for use in fully developed turbulence, as specified above. However, 
these requirements may be met by the use of silicon strip detectors as optical imaging elements 
in a particle tracking system. The strip detectors used in our experiment were developed to mea- 
sure sub-atomic particle tracks in the vertex detector of the CLEO III experiment operating at 
the Cornell Electron-Positron Colhder (Skubic et al. 1998). When applied to particle tracking 
in turbulence each detector measures a one dimensional projection of the image of the tracer 
particles. Using a data acquisition system designed specifically for the turbulence experiment 
several detectors may be simultaneously read out at 70,000 frames per second, making it possi- 
ble to measure two-dimensional or three-dimensional particle trajectories with very high spatial 
and time resolution (Voth et al. 2001; La Porta et al. 2001). 

The technical demands of particle tracking in turbulence and the scope of phenomena which 
become accessible when these demands are met can best be appreciated by examining one of the 
more dramatic particle trajectories recorded with the strip detector particle tracking system. The 
three-dimensional trajectory shown in figure 1 was recorded in a flow between counter-rotating 
disks at R\ = 970 (described below in Section 3) having an rms velocity of approximately 
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Figure 1 . Trajectory of a 46 /im diameter tracer particle in turbulence at Rx ~ 970 recorded at frame rate 
70,000 fps. The position of the particle at each of 278 frames is represented as a sphere. The acceleration 
magnitude is represented as the colour of the trajectory, as indicated by the scale. 

1 ms~^, a Kolmogorov time t,, of approximately 1/3 ms, and a Kolmogorov distance scale 77 of 
20 nm. 

The trajectory was recorded at 70,000 frames per second using a pair of strip detectors collect- 
ing primary and conjugate charge, as will be described in Section 2 below. The particle enters 
the measurement volume near the top right of figure 1 and appears to be trapped in a vortical 
structure. In the third, tightest turn of the helical motion the acceleration of the particle rises 
to 16,000 ms^^ within a time interval of 0.5 ms 1.5t^). During this trajectory the velocity 
components of the particle, shown in figure 2, oscillate wildly, spanning a range considerably 
larger than the rms velocity within a Kolmogorov time. The extreme fluctuations in velocity and 
acceleration, which occur in time-scales of order of the Kolmogorov time, make great demands 
on the particle tracking system. 

The trajectory in figure 1 and 2 is also noteworthy in that it seems to lie outside the charac- 
teristic range of the energy cascade. For the first part of the trajectory the particle appears to 
be caught an a very intense vortical structure. The tightest loop in the helical trajectory of the 
particle is about 300 /im in diameter, or about 15 ry, which is in the inertial subrange quite close 
to the dissipative scale of the turbulence. The period of the motion appears to be less than 1 ms, 
of the order of the Kolmogorov time r^. Yet the velocity fluctuations observed in this trajectory 
exceed the rms velocity, which would normally be associated with the largest scales of the tur- 
bulence (the so called energy-containing range). Events of this nature, which up until now were 
experimentally inaccessible, may be interpreted as a manifestation of turbulent intermittency. 

The remainder of the paper is organized as follows. In Section 2 the operating principles and 
capabilities of the silicon strip detector based particle tracking system are discussed. This in- 
cludes discussion of the optical system used to image the tracer particles on the strip detector 
(Section 2.2). In Section 3 the turbulent water flow between counter-rotating disks is described, 
and characterized using the silicon strip detector This entails measurement of the velocity statis- 
tics and estimate of the energy dissipation from the transverse velocity structure function. 

The main results are the investigation of the statistics of particle accelerations in turbulence for 
140 < R\ < 970, given in Section 4. (A brief account of part of this research has previously been 
published in Nature (La Porta et al. 2001).) In Section 4. 1 the probability density function (PDF) 
of the acceleration component is measured as a function of Reynolds number and component 
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Figure 2. Three components of the velocity of the particle shown in figure 2. 



direction. In Section 4.2 the scaling of the acceleration component variance is compared with 
the Heisenberg-Yaglom prediction. It is found that data is consistent with DNS results at low 
Reynolds number, and with the Heisenberg-Yaglom predicted scaling for R\ > 500. The flatness 
of the acceleration component is found to be quite large and to increase with Reynolds number, 
exceeding 60 at R\ = 970. The acceleration component autocorrelation function is shown in 
Section 4.4. The autocorrelation function is shown to scale with the Kolmogorov time, and to 
cross zero at a time of approximately 2.1 r^. Finally, the influence of particle size and density are 
studied in Section 5. It is found that the acceleration of relatively large particles (diameter > 5 r^) 
can be significantly smaller than that of infinitesimal fluid particles, and that, in this experiment, 
particle size and not density is primarily responsible for this. However, the 46 /im particles used 
for the acceleration measurements differ from ideal fluid particles by only a few percent. Two 
appendices are also included. The first details the algorithms used for extracting particle tracks 
from the strip detector output, and the second compares the results published here with previous 
measurements, which were made with a conventional detection system which had much lower 
resolution. 



2. Particle Tracking System 

2.1. Strip Detector 

A brief description of the the strip detector tracking system is provided below. It is described in 
detail in (Voth et al. 2001). 

The silicon strip detector is essentially a large planar photodiode (sensitive area 5 1 mm x 
25.6 mm) which is segmented into 512 sense strips, as shown in figure 3. Optical radiation in- 
cident on the detector creates electron-hole pairs, and the holes are collected by p-type strips 
patterned on the front surface of the n-type detector. Charge is conducted off the detector chip 
through metallic leads evaporated onto the surface of the detector which are oriented perpen- 
dicular to the sense strips. The positive charge collected by the array of strips gives a one- 
dimensional projection of the light intensity incident on the detector, allowing the x coordinate 
of a particle that is imaged onto the detector to be measured. However, for certain settings of 
the detector bias, electrons appear to become trapped in meta-stable surface states and couple 
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Figure 3. Simplified scliematic of the strip detector. Dark grey bars represent p-type sense strips. Liglit 
grey bars superimposed on tlie darlc bars represent metal leads evaporated on to the surface of the detector. 
The small circles along the diagonal represent paths that connect sense strips to the leads. Under normal 
circumstances, only "primary" charge is collected indicating the x coordinate of light spots incident on the 
detector. However, the detector bias may be adjusted so that "conjugate" charge is also detected by the metal 
leads indicating the z coordinate, as described in the text. 

capacitively to the metallic leads. Under these circumstances, negative "conjugate" charge in- 
dicates the projection of the intensity on the axis defined by the leads (which is perpendicular 
to the axis defined by the sense strips). By distinguishing the primary (positive) and conjugate 
(negative) peaks the two-dimensional projection of a particle trajectory may be measured using a 
single detector. Positions measured from conjugate charge have greater uncertainty, so for most 
of the data presented below, and in particular for the acceleration measurements, primary charge 
was used for position measurement. 

The 512 strips of each strip detector are connected to an array of four resistor-capacitor (RC) 
chips which provide an individual bias resistor and coupling capacitor for each channel. The 
output capacitors are connected to an array of four integrated amplifier/multiplexer chips, which 
provide a shaper amplifier and sample-and-hold amplifier for each strip. After a frame of data is 
latched on the array of sample and hold amplifiers, analog multiplexers integrated into the am- 
plifier chips output the strip intensities as four 10 Msample/sec waveforms. A detector controller 
generates frame-readout triggers for the strip detectors, readout electronics and an acousto-optic 
modulator, and may be configured for readout rates in the range 5 kHz-70 kHz. The acousto-optic 
modulator is used to strobe the illumination before the readout of each frame, which is neces- 
sary to optimally drive the shaper amplifiers. The output waveforms are captured using a pair of 
dual-channel digital oscilloscope boards mounted in a PCI bus slave computer The oscilloscope 
boards can store 4,000 frames of data in internal memory, which is subsequently downloaded 
into the computer's main memory. The slave computer performs pedestal subtraction, thresh- 
olding, and stores compressed data to a local hard disk concurrently with the acquisition of the 
next frame of data. The maximum duration of continuous acquisition ranges from 800 ms at 
5 kHz to 58 ms at 70 kHz. The time required to process and store one 4000 frame sequence is 
approximately 1 s, and this limits the duty cycle of the acquisition system at high frame rates. 

2.2. Optical Imaging System 

The optical configuration used to image particles is shown in figure 5. The illumination beam, 
generated by a 6 W continuous wave Argon-Ion laser and gated by the acousto-optic modulator, 
is directed through a glass window and passes through the centre of the flow chamber. The beam 
is a TMqo Gaussian mode and a beam expander is used to obtain a spot radius oj k, 1.0mm with 
negligible divergence. This beam radius is chosen to so that the beam fills the field of view of 
the detector. Acceleration measurements were taken in the configuration shown in figure 5(a), 
in which two strip detectors view a common image and measure two-dimensional coordinates. 
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Figure 4. The readout cluster. Each strip detector is read out using an individual slave computer A field 
programmable gate array (FPGA) detector controller generates timing signals for the strip detector data 
acquisition hardware, and acousto-optic modulator. A master computer controls the slave computers via 
Ethernet messaging and provides the user interface for the system. 



Some three-dimensional tracks were recorded in the configuration shown in figure 5(b), in which 
two detectors view the volume from different view-ports so that their detection volumes overlap. 

The design of the image system, sketched in figure 5 is determined by the requirement that a 
small measurement volume be imaged at high magnification by optics which are placed outside 
the turbulence chamber These design requirements were met by the two stage imaging system 
shown in figure 6(a). In order to maintain acceptable depth of field, the aperture shown after the 
first imaging lens (LI in figure 6) is used to restrict the numerical aperture of the imaging system 
to a value of 0.03, giving a spot size w 3G0/im for particles throughout the measurement volume. 
As a result, at least three strips were illuminated by each particle image, making it possible to 
locate the particle with sub-strip resolution by fitting each peak to a Gaussian function. 

Although the magnification of the optical system defined in figure 6 is nominally f^/ f2, it is 
useful to vary the configuration of the system in order to adjust the magnification and position 
of the focal plane. In particular, Li may be repositioned to change the magnification of the first 
stage of the system, and is replaced by a pair of lenses whose separation dependent effective 
focal length replaces /a. For the acceleration data presented below, Li is a compensated doublet 
with focal length 15 cm, L2 is a doublet of focal length 3.8 cm and L3 is replaced by a pair of 
lenses consisting of a 30 cm focal length piano convex lens and a —50 cm piano concave lens 
separated by approximately 10 cm. The detector is placed approximately 25 cm behind the piano 
concave lens. The magnification of the system is approximately 12.8, so that the 100 /im pitch 
of the strip detector corresponds to approximately 7.8 /im in the fluid volume, giving a field of 
view of 2.00 mm x4.00 mm. The effective magnification of the system is measured to within 2% 
using a calibration target that is positioned in the active region at the centre of the flow chamber. 
(The variation of the magnification is less than 1 % over the illuminated volume, so it was not 
necessary to express the magnification as a function of depth.) 

For acceleration data, the strip detectors are biased such that conjugate charge is suppressed 
allowing them to measure a single coordinate using the primary charge signal. In order to measure 
two-dimensional trajectories, a beam-splitter is used to project the same image onto two detectors 
(see figure 5(a)), which are oriented to measure orthogonal coordinates. Because of the 2: 1 aspect 
ratio of the detector, the fields of view of the two detectors do not coincide, as shown in figure 7. 
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Figure 5. (a) Optical layout for acceleration measurements, viewed along the axis of the cylindrical tur- 
bulence chamber (z). The illumination beam is gated by an acousto-optic modulation (AOM) before being 
passed through the centre of the turbulence chamber Optics are used to image the central volume through a 
view-port at 45° with respect to the illumination. A beam-splitter allows the image to be projected on two 
strip detectors which are oriented to measure two orthogonal coordinates. An additional view-port may also 
be used to measure the third coordinate, (b) Optical layout for three-dimensional particle measurement. The 
apparatus is the same except that imaging optics are present on both ports. In this case strip detectors are 
operated with conjugate peaks enabled (see text) so each port has a detector which measures two position 
coordinates, on each port. 

In this case, the full area of each detector can be used to measure one-dimensional trajectories, 
but two-dimensional trajectories can only be measured in the region where the two fields overlap. 
As a result, only half of the 512 strips are available when analyzing two-dimensional trajectories. 

2.3. Determination and Characterization of Tracks 

Data files stored by the experiment consist of thresholded intensity data simultaneously acquired 
by the two strip detectors. These data files are passed through several stages of analysis. This 
task is similar in principal to that employed in standard particle tracking studies using CCD 
cameras, but the algorithms employed must be tailored to the pecuUar characteristics of one- 
dimensional projection images. A description of the algorithm used for track extraction is given 
in Appendix A. Secondary data files, consisting of lists of tracks measured for each data sequence 
are stored for subsequent analysis. 

The accuracy of the particle positions is estimated by fitting a straight line to very short track 
segments (shorter than -r,,) and measuring the mean deviation between the data points and the fit. 
The deviation is found to depend most significantly on the peak height. Figure 9 shows the error 
as a function of peak height for a typical run. For extremely weak peaks (maximum intensity < 
0.02 V) the error can be of order 1/2 strip, but the the analysis routines are configured to limit the 
data to larger intensity, so that the mean error is < 0.1 strips. The error for individual trajectories 
can vary depending on focus and aperture, but always remains below 0.2 strips. 

One can also consider the mean deviation of the peak from the linear fit as a function of the 
position (in strips) modulo 1, shown in figure 10. The dependence of the mean deviation on the 
position indicates a nonlinearity in the interpolation of the peak centres which is smaller than, but 
of the same order as the random uncertainty. This nonlinearity is compensated for by subtraction 
of the mean deviation from the peak positions before processing of the trajectories. 

The kinematic properties of the tracks are measured by polynomial fits (parabolic for acceler- 
ation, linear for velocity) to the position vs time data. The fits are made using the standard least 
squares algorithm, with the relative weight of each data point proportional to the inverse square 
of the estimated error (as a function of peak intensity, figure 9). The length of track over which 
the fit is performed is a complex issue, which is discussed in detail below. 
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Figure 6. A simplified schematic of the optical system is shown. Plano-convex lens Li is place just outside 
the fluid chamber and is used to image the active volume with 1:1 magnification. A short focal length lens 
L2 is positioned so that its front focal plane (/2) lies at the centre of the image formed by Li. As a result, 
point sources in the image are transformed into parallel ray bundles by L2. The bundles are focused by 
L3 at its rear focal plane (/s). A beam splitter may be inserted between L2 and L^, to allow the image 
to be projected on a second detector, as shown in figure 5. The magnification of the system is given by 
the ratio of the focal lengths fs/ f2 - The light collection of the imaging system is controlled by an aperture 
placed behind lens Li . (b) The tracer particles are transparent polystyrene spheres with density 1 .06 g cm~"^ 
and diameter 46 /im and are detected via specular reflection from the internal and external surfaces. With 
illumination at 45° and polarization perpendicular to the plane of incidence, the internal reflection is much 
stronger than the external reflection so that the particle appears as a single point-like source of light. 

One of the greatest challenges in the measurement of accurate statistics of Lagrangian vari- 
ables is the control of sample biases. From a physical standpoint, one must take care that the 
tracer particles uniformly sample the fluid volume and follow the flow field. This requires that 
the particles be sufficiently small and that the density match between the particle and the fluid 
be sufficiently close. This issue is addressed in Section 5, below. It is also necessary to ensure 
that the measurement and analysis procedures do not introduce biases in the measurements of 
the particle trajectories. For instance, if one were to make one measurement for each particle that 
enters the measurement volume, a distorted velocity distribution would be obtained, since the 
rate at which particles enter the measurement volume is itself proportional to the velocity of the 
particle (Buchave et al. 1979; Voth et al. 1998). However, the time that a particle will remain in 
the measurement volume is inversely proportional to the velocity, and these two factors cancel. 
Ideally, one would achieve uniform sampling within the measurement volume by continuously 
measuring the kinematic properties of the track from the time that a particle enters the volume 
to the time that it exits. In practice this is impossible to achieve because the acceleration and 
velocity are measured by fitting to a polynomial function, and so the variables can not be mea- 
sured until the particle has been in the measurement volume for a finite time. Another difficulty 
is that measurement is not possible when particles cross paths or traverse inoperative pixels. The 
strategy employed is to measure the variables as many times as possible along the trajectory, 
and make the total statistical weight of these measurements proportional to the total length of 
the track. This seems to give the best approximation of uniform sampling of the measurement 
volume. 

3. Characterization of Flow 

The goal of this study is to explore universal characteristics of turbulent flows. However, our 
flow deviates significantly from the ideal of homogeneity and isotropy. The apparatus used to 
generate the turbulent water flow is described in Section 3.1. Standard techniques for character- 
izing turbulence using anemometers are not feasible. Therefore, we have used the strip detector 
particle tracking system to characterize the flow in terms of standard Eulerian quantities. This 
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Figure 7. The relative fields of view of two strip detectors oriented as in figure 5(a) to measure perpen- 
dicular coordinates, where the grey stripes represent the orientation of the charge collecting sense strips. 
Because of the 2: 1 aspect ratio of the detector, only half of the field of each detector overlaps with the field 
of the other, so that only 256 strips can be used for coincident measurements. 
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Figure 8. (a) Raw data from the two strip detectors. Only the 256 overlapping strip segments are shown. 

(b) Positions vs time measured from the raw data, where gaps are due to the presence of inoperative strips. 

(c) Reconstruction of the two-dimensional track from the position data in (b), where gaps have been bridged 
by linear interpolation. 
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Figure 9. Estimate of particle position error. The graph shows tiie rms deviation of tiie peak centre from a 

linear fit to a short segment of ttack. 



includes measurement of the scaling of the rms velocity components with propeller speed (Sec- 
tion 3.2), an estimate of the rate of energy dissipation in terms of the velocity structure functions 
(Section 3.3), and a measurement of the anisotropy of the mean velocity field at the centre of the 
flow (Section 3.4). 

3.1. Turbulence Generator 

The turbulence is generated in a flow of water between counter-rotating disks in a cylindrical 
container (Voth et at. 1998; La Porta et at. 2000), as represented schematically in figure 1 1 and 
shown in figure 12. The container, mounted vertically, is 48.3 cm in diameter and 60.5 cm long 
with 8 planar windows mounted flush to the surface of the cylinder at equal angles along ihe mid- 
line of the chamber. All physical quantities are defined with reference to a Cartesian coordinate 
system in which the z axis corresponds to the axis of symmetry of the cylinder and the x and 
y axes correspond to the optical axes of the two imaging systems represented in figure 5(b). 
By symmetry the two transverse coordinates (x and y) are equivalent, and distinct from the 
axial coordinate (z). The disk shaped propellers are open-ended cylinders 20 cm in diameter and 
4.3 cm deep with twelve internally mounted radial vanes. The propellers are spaced 33 cm apart 
and are each driven by 0.9 kW computer controlled controlled dc motors that are coupled to the 
propellers with variable speed reducers. A smaller cylindrical tube surrounds each propeller and 
stationary radial vanes between this inner cylinder and the container wall have been installed to 
inhibit large scale rotation of the flow. For the studies described below, the propeller rotation rate 
is varied from 0.15 Hz to 7.0 Hz. The lower Umit is set by the motors and speed reducers and 
the higher limit is set by the temporal resolution of the detector. (The detectors are not able to 
adequately resolve the Kolmogorov time at the maximum propeller speed of 9 Hz.) 

The averaged flow produced by the propellers can be interpreted as a superposition of two 
basic components, a pumping mode and a shearing mode. Centrifugal pumping by the propellers 
produces the flow represented schematically in figure 11(a). The resulting mean strain field at 
the centre of the chamber (represented by the arrows in figure 1 1(a)) tends to enhance the axial 
component of the vorticity. In addition, fluid near the top and bottom of the cylinder tends to 
rotate collectively with the counter-rotating propellers, creating a shear layer around the edge of 
the flow midway between the propellers, as represented in figure 1 1(b). Particle accelerations are 
measured in a 4 mm^ volume at the centre of the flow chamber. Quantitative measures of the 



Measurement of Particle Accelerations in Fully Developed Turbulence 



11 



^0.10 




0.0 

position (modulo 1) 

Figure 10. Deviation of particle position from fit as a function of position modulo 1. The circle indicates 
the mean and the error bars indicate the standard deviation of the deviation of the peaks from the fit. 



flow, described in detail below, are derived from measurements in this volume, and additional 
measurements in a larger 15 mm x 30 mm volume. More details regarding this flow are given in 
references. (Voth et al. 1998) and (La Porta et al. 2000). (See also Appendix B.) 

3.2. Measurement of the Velocity Fluctuations 

Velocity statistics may be obtained from analysis of matched trajectories obtained with the appa- 
ratus in its standard configuration, shown in figure 5(a). As mentioned above, in order to obtain a 
correct estimate of the rms velocity, it is necessary to take care to sample trajectories uniformly 
by insuring they are continuously sampled as long as they remain in the measurement volume. 
The PDF of the velocity is shown for linear and log scales in figure 13. It is evident that the stan- 
dard deviation of the transverse component exceeds that of the axial component by about 50%. 
The distributions for both components are approximately Gaussian; the flatness is 2.8 for the ax- 
ial and 3.2 for the longitudinal component. These flatness values are independent of the propeller 
rotation frequency (NouUez et al. 1997). It may be noted that there is a small dip near zero veloc- 
ity in the axial velocity component PDF. This occurs because the light sensitivity of the detector 
decreases when a peak remains on the same pixel for consecutive frames. (This is due to an inef- 
ficiency in the shaping amplifiers used for this detector.) As a result, the measurement volume is 
effectively smaller for trajectories with near-zero velocity, causing a measurement bias. It should 
be noted that measurement biases depend mostly on the velocity and do not affect acceleration 
measurements, except to the extent that the acceleration and velocity are correlated. 

The scaling of the rms velocity with propeller rotation frequency is shown in figure 14. The ex- 
pected linear dependence on the frequency is observed for both components, which is consistent 
with the assumption that the nature of the large scale flow is independent of the stirring velocity 
over the range of Reynolds number studied. The inset shows that the deviation from the linear 
scaling law does not exceed a few percent. This deviation is most pronounced at low propeller 
speeds, perhaps indicating that the turbulence is not "fully developed" at the lower end of the 
Reynolds number range. 

The ratio of the rms transverse velocity to the rms axial velocity is shown in figure 15. The 
ratio varies only a few percent over the full range of propeller speeds, again indicating that the 
large scale structure of the flow does not change as the Reynolds number is varied. 
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Figure 1 1 . Schematic representation of the flow between counter-rotating disks decomposed into (a) the 

pumping mode and (b) the shearing mode. 



3.3. Energy Dissipation 

In order to compare results with Kolmogorov scaling predictions, it is essential to measure the 
energy dissipation rate e, since all statistical quantities are assumed to depend on this quantity 
and upon the kinematic viscosity v. The energy dissipation e is given by 



e = 2v{sijSij) 

where is the fluctuating rate of strain tensor, defined by 



duj 



duj 
dxi 



(3.1) 



(3.2) 



where m = Ui — {Ui) is a component of the velocity fluctuation (Pope 2000). Unfortunately 
Lagrangian particle tracking does not allow us to measure velocity gradients, so a direct mea- 
surement of the energy dissipation is not possible. 

There are several possible indirect methods for measuring the energy dissipation from particle 
tracking data (Ott & Mann 2000). The most accessible method requires measurement of second 
(or thkd) order velocity structure functions. The velocity structure functions are the moments of 
the velocity differences between two points separated by a fixed difference. The longitudinal and 
transverse structure functions are calculated from velocity components parallel to and perpendic- 
ular to the line separating the two points, respectively. From the second order structure functions, 
the energy dissipation can be obtained by comparison with the K41 scaUng relations, 

C2(er)2/3 



D 



LL 



(3.3) 
(3.4) 



where L and TV designate longitudinal and transverse, respectively, and the separation r is as- 
sumed to be within the inertial subrange. C2 is an approximately universal constant that has been 
determined empirically (Monin & Yaglom 1975; Pope 2000). 

All of these structure functions can be calculated from three-dimensional Lagrangian trajec- 
tory data by simultaneously measuring the velocities of pairs of particles. However, the particle 
tracking system used in this experiment is currently not capable of three-dimensional tracking in 
a volume which would encompass inertial range particle separations due to power limitations of 
the Argon-Ion laser used for illumination. To overcome this, we have developed a technique for 
measuring the transverse second order velocity structure function using two-dimensional particle 
tracking in conjunction with a light sheet. 
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Figure 12. Apparatus, consisting of turbulence chamber with counter-rotating disks (foreground), 
illumination beam entering from foreground and two detectors (background). 

The structure function measurements are performed using the configuration shown in fig- 
ure 16. There are several significant differences between this configuration and the standard 
configuration used for single particle acceleration and velocity measurements. The magnifica- 
tion has been reduced to 2.89, giving a field of view of 17.7 mm (> SOOry), which allows inertial 
range particle separations to be observed. The optical system is schematically similar to figure 5, 
except that the beam-splitter is omitted and a single detector is used to measure tracks. The bias 
of the detector is set so that conjugate peaks are enabled and two-dimensional trajectories can be 
obtained from the single detector output, as discussed in Section 2.1 above. In order to maintain 
adequate illumination intensity over such a large field of view, the illumination beam is config- 
ured as a light sheet approximately 0.1 mm thick and 10 mm wide. The Ught sheet was created 
using a standard Galilean telescope to expand the beam to a large TMqo mode, then a second 
cylindrical telescope to compress the horizontal axis and create an elliptical mode. The measure- 
ment would be simpler in principle if the Ught sheet were parallel to the image plane, but in prac- 
tice it must be oriented at 45° with respect to the optical axis in order to obtain sufficiently strong 
light scattering from the particles. The data obtained from this configuration consists mostly of 
short tracks which are created as the particles pass through the light sheet. The depth of the light 
sheet is chosen so that tracks have sufficient length for an accurate velocity measurement to be 
made. 

Velocity structure functions are measured as a function of r by calculating velocity differences 
for all coincident pairs of particles and compiling statistics under the condition that the separation 
distance hes within an adjustable range. Such a pair of particles will appear to the detector as 
shown in figure 16(b), and it is possible to measure two velocity components for each particle 
and two components of the particle separation vector. However, as illustrated in figure 16(a), the 
separation vector does not he in the same plane as the velocity components. Using the coordinate 
system defined in figure 16 the coordinates x and z are measured explicitly and we may assume 
that y = X. The velocities are determined by tracking the particle for the short time that it 
remains within the finite thickness of the light sheet, and gives the projection of the velocity on 
the imaging plane. We therefore measure the components Ux and Uz, but Uy is unknown. 
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Figure 1 3 . (a) The PDF for the velocity, (b) Same data as (a) plotted on a semi-log scale. 



1.0 



0.5 

3 



0.0 









T 1 


1 ' 


1 1— 


- 1 .00 




i 


















■ 0.98 


o 










■ 0.96 


?^ 




o transvers? 










& axial 






' 0.94 


A 










( 




2 


4 6 

















transverse 










A 


axial 



2 4 6 

propeller speed (Hz) 



Figure 14. The axial and transverse components of the rms velocity as a function of propeller speed. The 
lines indicate the best fit to a linear relationship Wj = fc,/, with ka = 0.0910 for the axial component and 
kt — 0.137 for the transverse component. The inset shows the rms velocity components normalized by the 
linear scaling law. 
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Figure 15. The ratio of transverse (x) to axial (z) components of the rms velocity. 
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(a) view of experiment from above 


(b) detector view 
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Figure 16. (a) Experimental setup for dissipation measurements viewed from above (along the axial direc- 
tion) showing the orientation of the 45° light sheet which illuminates the detection volume. (b)Illustration of 
detector view for dissipation measurements. The detector sees the particle separation and velocity projected 
on the X — z plane, where z is the axial coordinate. 



Using this geometry, it is not possible to measure the longitudinal structure function because it 
is in general not possible to measure a velocity component along the separation vector. However, 
it is always possible to measure one of the velocity components perpendicular to the separation 
vector. To find the measurable transverse velocity component, we make use of the fact that this 
component must be perpendicular to both the separation vector Vyi and to the vector normal to 
the image plane (y). Such a vector may be constructed by taking the cross product of ri2 and y. 
The transverse velocity component is then 

M_L2 - = (U2 - ui) • (ri2 X 2/)/||ri2 X yll, (3.5) 

where Ui and U2 are the velocities of the two particles. The dot product can be evaluated even 
though the y component is unknown because the y component of Vyi x y is identically zero. 

The structure functions are calculated from the second moment of (u_l2 — wj^i) conditional on 
ri2, 

i?NN(r) = ((^^±2 - V < ||ri2|| < (1 + s)r^ (3.6) 

where the relative bin width s is maintained constant as r is varied and where f is the mean value 
of ||ri2|| for all events which satisfy the condition. 

Raw data used in the structure function are shown in figure 17(a) and x and z coordinate 
tracks (calculated from primary and conjugate peaks, respectively) are shown in figure 17(b). The 
matching of x and z tracks is performed using intensity correlations (as described in Appendix A 
below) The results of the matching are shown in figure 17(c). The scaling of -Dnn with r is 
shown in figure 18. Both exhibit the expected r^/^^ scaling over a substantial range. The scaling 
range is limited at small r because particles are too close together to be well localized by the light 
sheet and at large r because their separation is approaching the integral length scale. The energy 
dissipation rates obtained from the fits in figure 18 are shown in Table 1, and and are consistent 
with the expected e = u^/L scaling with L = 7.1 cm. Once the energy dissipation has been 
determined, it is possible to calculate the Taylor microscale, A — {Ibvv? /eY/"^ and the Taylor 
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Figure 17. (a) Raw intensity data, in which dark areas indicate primary (positive) charge and light areas 
indicate conjugate (negative) charge, (b) Position vs time for x (primary) and z (conjugate) coordinates, (c) 
Reconstruction of x-z trajectories after track matching. The dashed and dotted line circles indicate a pair of 
simultaneous tracks from which velocity differences may be measured. 



microscale Reynolds number 

R\ = 



V 



^5l/2gl/6j^2/3 



(3.7) 



3.4. Mean Rate of Strain 

The optical configuration used for the dissipation measurements is also useful for investigation of 
the structure of the large scale flow near the centre of the apparatus. Due to the symmetry at the 
centre of the flow, the off-diagonal elements of the mean rate of strain tensor §^ are zero, and 
the diagonal elements are constrained by symmetry in the transverse plane and the incompress- 
ibility condition, = (using the summation convention), so that = —2-^^ = —2-^^. 

Measurements at 5.0 Hz and 2.5 Hz show a linear relationship — m/O. 0492m. This may 
be compared with a component of the fluctuating strain, which is related to the dissipation in 

an isotropic flow by e 
therefore 

m (i5.L)v^ 



15u (^f^)^- The ratio of the fluctuating strain to the mean strain is 



K dz 



0.0492m 



(3.8) 



This implies that at the maximum propeller speed of 7 Hz, the mean strain is about 2% of the 
fluctuating strain and rises to about 15% of the fluctuating strain at the minimum propeller speed 
of 0.15 Hz. 



4. Results: Particle Acceleration Measurements 

The primary subject of this paper is the study of fluid particle accelerations in fully devel- 
oped turbulence. In contrast to the fluid particle velocity, which is the same quantity that would 
be measured by a fixed probe at the same location, the particle acceleration can only be mea- 
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Figure 18. Log plot of the scaling of Dnn with r. The square symbols indicate data at 5 Hz and the 
round symbols indicate data at 2.5 Hz. The straight lines indicate best fit to as indicated in Table 1 . 



sured using Lagrangian techniques. The acceleration of a fluid particle a+ corresponds to the 
substantive derivative of the velocity 

a+^^ + (u.V)u. (4.1) 

In order to determine a+ from Eulerian measurements, it would be necessary to know du/dt as 
well u and Vu at a point in space (which is possible in DNS but not in experiments). In terms of 
the fluid particle acceleration, the Navier-Stokes equation is 

, VP . 

a+ = + iyV^u, (4.2) 

P 

where P is the pressure and p is the fluid density. In fully developed turbulence, the viscous term 
is small compared with the pressure gradient term, so a measurement of a+ gives information 
about the pressure gradient, which is difficult to measure experimentally. 

Measurements of the acceleration require the lowest possible position error, and were made ex- 
clusively from primary charge readout (rather than from less accurate conjugate charge readout). 
Accelerations were measured from two-dimensional trajectories recorded in the configuration 
illustrated in figure 5(a). Accelerations are therefore known in the x-z plane, giving one axial 
and one transverse component. (Due to symmetry, the statistical properties of the unmeasured y 
coordinate are expected to be identical to those of the x coordinate.) 

4. 1 . Acceleration Distribution 

The PDF of the x (transverse) component of the acceleration is shown on linear and logarithmic 
scales for several values of R\ in figure 19. It is found that the distributions have a stretched 
exponential form for all measured values of Rx, but that the extension of the tails increases with 
R.\. The distributions are plotted for fits to a finite time interval r/ of 0.75t,j. Although there is 
no qualitative change in the distribution as the fit time is varied over a range 0.5r,, < t/ < 2t,,, 
the moments of the distribution tend to increase as the tj is reduced. Estimates of the variance 
and flatness of the acceleration distribution ((a^) and {af)/{af)^, respectively) must therefore 
be obtained by measuring these quantities as a function of ry and extrapolating to zero, as will 
be described in Section 4.2 below. However, on the basis of figure 19 it is obvious that the tails 
extend far beyond those of a Gaussian distribution of the same variance and that the flatness of the 
acceleration is very large (> 40, compared with 3 for a Gaussian). This result is consistent with 
the large pressure gradient flatness values measured at low Reynolds numbers in DNS(Vedula & 
Yeung 1999). 
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f u e L 

(Hz) (ms-^) (m^s-^) (m) 



2.5 0.3095 0.406 0.073 
5.0 0.6190 3.367 0.070 

Table 1 . Turbulence parameters for dissipation data. The integral length scale L is calculated using 
e = w^/L and C2 = 2.13(Sreemvasan 1995). 



The acceleration component PDF may be parameterized by the phenomenological function 

P(a) = Cexp(-[a2]/[(l + ja/J/aDa^]), (4.3) 

where /3 = 0.539, 7 = 1.588, a = 0.508 and C = 0.786 was obtained for the data at Rx = 970. 
It has been shown by Holzer & Siggia ( 1 993) that if the tails of the velocity difference distribution 
are exponential, then the acceleration PDF should exhibit scaling P{a) oc exp(— a^/^). This may 
be compared with the tails of equation 4.3 which have the form exp( — |a|°'^^). 

As a result of the long tails of the acceleration PDF, very large amounts of data are required for 
convergence of the fourth moment, which is required for the calculation of the flatness. Figure 20 
shows the contribution to the second and fourth moments as a function of acceleration. It is 
apparent that even with 5x10^ acceleration measurements, the convergence of the fourth moment 
is marginal. 

The extent to which the anisotropy of the flow affects the acceleration PDF is illustrated in 
figure 21. The plot shows the PDF's for and at Rx — 970. Comparison with figure 13 
shows that the acceleration is much less anisottopic than the velocity. The small difference in the 
variances of Ux and wiU be discussed below in Section 4.2. 

4.2. Acceleration Variance 

The variance of particle accelerations in a turbulent flow was first predicted on the basis of the 
1941 scaling theory of Kolmogorov (Kolmogorov 1941b, a) by Heisenberg and Yaglom (Heisen- 
berg 1948; Yaglom 1949). The variance of the acceleration components is given by 

{aiaj) = aoe^^^p-'/^Sij (4.4) 

where oq is predicted to be a universal constant which is approximately 1 in a model assuming 
Gaussian fluctuations (Heisenberg 1948). The form of this scaling law can be deduced from the 
assumption that the acceleration is a dissipation scale quantity, and must be determined only by 
e and i^. 

Deviations from the Heisenberg- Yaglom scaling law are expected to arise from turbulent in- 
termittency. Using the refined similarity theory, equation 4.4 is replaced with 

{a,a,) - aou-''^^6ij (e^/^^ , (4.5) 

where e^, the energy dissipation averaged over a sphere of radius r, has taken the place of the 
mean energy dissipation, e. Using the log-normal model for the moments of e^, this yields 

{aiaj) oc {L/nf^'^ oc rI^''^^ = R^^^ (4.6) 

where a value of 1/4 has been used for the intermittency exponent, 112. Other models of intermit- 
tency have been developed, such as the explicitly Lagrangian model of Borgas (Borgas 1993), 
which predicts ao oc Rx^'^^^- 
Direct numerical simulation of turbulence has shown that at low Reynolds number ao Rx^^^ 
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Table 2. Turbulence parameters for acceleration data with v = 9.89 x 10 ''vo? s ^ (water at 20.6°). / is 
the propeller rotation frequency and At is the strip detector frame period, u — ((u^ + + u^^)/?))^^'^ is 
the rms velocity. The sweeping time is calculated from u and the 2.00 mm field of view of the detector. 
The energy dissipation is calculated from e = /L with L = 0.071m, as discussed in Section 3.3. The 
Taylor microscale Reynolds number Rx is calculated from equation 3.7 and the classical Reynolds number 
Re is defined by Re = Ra^IS. 
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Figure 19. The PDF of the transverse (x) component of the acceleration at three values of R\ plotted on 
(a) linear and (b) log scale. Each acceleration distribution is measured by parabolic fits over 0.75t^ and is 
normalized by its standard deviation. The dashed curve is a Gaussian with the same standard deviation and 
the solid curve is the parameterization defined in equation 4.3. 
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Figure 20. The curves show the relative contribution to the second and fourth moments of the transverse 
acceleration component, (o^P(aa;)/(a^) and {a'iP{ax) / ifli)) respectively, as a function of acceleration. 
P{ai) is the PDF of the i component. 



(or e^/^^, using equation 3.8 and assuming constant L and v). This is equivalent to an overall 
scaling of the acceleration variance with e^^/^^, which is a relatively small deviation from the 
Heisenberg-Yaglom prediction of e^/^. 

In principal, the acceleration variance may be calculated by taking the second moment of the 
distributions shown in figure 19. However, in order to calculate the particle accelerations the 
track must be fit over a finite time, and the moments of the distribution depend on this fit time 
interval. The issue is illustrated in figure 22, which shows a typical particle trajectory. The raw 
trajectory appears to be straight, but when the mean velocity is subtracted off, the particle is seen 
to be undergoing a small time varying acceleration, which causes it to deviate from a straight 
trajectory by a distance corresponding to a few strips on the detector. The acceleration of this 
track is approximately 90 ms~^, which is about one quarter of the rms acceleration — a typical 
value. Clearly the correct acceleration will be measured from this track only if a parabolic fit is 
made over an appropriate time interval. The fit over 150 frames does not conform to the trajectory 
and underestimates the acceleration. The fit over 5 frames conforms to the position measurement 
errors and dramatically overestimates the acceleration. The fit over 40 frames (1.8 r^) appears to 
conform to this particular trajectory. 

Ideally, we would like to fix the fit interval at a value where the noise is adequately averaged 
and yet the parabolic fits are able to conform to all of the particle trajectories. Figure 23 shows 
the normalized variance of the acceleration distribution as a function of the fit interval t/, and 
demonstrates that no value of r/ exists which satisfies this criteria, since there is no range of r/ 
where the acceleration variance is independent of tj. For r,, < tj < 9t,j there is an approxi- 
mately exponential dependence of the acceleration variance on the t/ which is due to the failure 
of the fits to conform to the true particle trajectories. For r/ < the acceleration variance rises 
dramatically with a w t^'' power law dependence. This is the t/ dependence which would be 
obtained from uncorrelated Gaussian distributed noise, and evidently arises from the position 
measurement error. 

The fact that the onset of position uncertainty occurs at a value of t/ where the fits fail to fully 
conform to the turbulent trajectories indicates that the frequency spectra of these two processes 
are not distinct, but overlap. There is therefore no way to distinguish these two contributions on 
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Figure 21 . Acceleration distributions for transverse (x) and axial {z) components of the particle accelera- 
tion for a run at Rx = 970. The acceleration was calculated for fits over 1 r,,. The standard deviations for 
X and z are 382 m and are 364 m respectively. 

any given track. However, we can separate the two effects by making use of the fact that the two 
contributions to the acceleration variance have very different scaling with tj. The procedure we 
use is to fit the measured uq (the normalized acceleration variance) to the function 

/(r) = At^ + Cexp(DT + Et^), (4.7) 

where A, B, C, D and E are fit parameters. The power law term represents the contribution 
from the position noise and the exponential term represents the contribution of the turbulence to 
the acceleration variance. (The term is added to model the slight deviation from exponential 
dependence observed at large r/.) The best estimate of ao is then obtained by evaluating the 
exponential term in the limit where r/ 0, so that ao = C. It is known that this extrapolation 
must overestimate the value of ao because the slope of the oq vs t/ curve must go to zero at 
Tf = since the tracks are differentiable. A simulation of the detection process, described in 
Section 4.5 below, indicates that this overestimate is 10% at R\ « 240, which is comparable to 
the random measurement error. All normalized acceleration variance data presented is rescaled 
to correct for this overestimation. This correction depends on the assumption that the dependence 
of the measured value of ao on (r/ /r^) has a universal form for small r/. 

The procedure described above was used to calculate the acceleration variance for x and z 
components of the acceleration for all values of the Reynolds number. Figures 24 (a) and (b) 
show ao as a function of r/ for the x and z components of the acceleration, respectively. The 
maximum value of rj for which the acceleration variance can be calculated is determined by the 
length of time the particles remain in view. This is estimated using the sweeping time Tg, which 
is the length of time a particle moving at the rms velocity would remain in the detection volume 
(tabulated in Table 3.4). For purposes of fitting to equation 4.7 the range of each ao curve in 
figure 24 was hmited at low r/ to the value where the power law term contributes an order of 
magnitude more than the exponential term, and at high t/ to 1.15ts. 

Certain conclusions may be drawn from figure 24. The curves for the three high Reynolds 
number runs collapse onto a single curve, indicating that the scahng of the acceleration with 
time is well approximated by K41 scahng. At low Reynolds numbers, the slopes of the curves 
and their extrapolations to rj = fail to collapse, indicating that K41 scaling has broken down. 
However, the interpretation is somewhat more complex because the measurement process brings 
an additional time scale into play. As the Reynolds number is reduced, the ratio of the residence 
time of the particles r, to the Kolmogorov time becomes smaller. If there is any correlation 
between residence time and acceleration, then this implies that similarity can not be achieved by 
rescaling the time axis with r^. By varying the size of the measurement volume (and hence the 
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Figure 22. Top panel shows an apparently straight particle trajectory recorded at Rx = 970. The duration 
of the track is approximately 150 frames, 2 ms, or 7 t,,. The bottom panel shows the same trajectory after 
the mean velocity has been subtracted off. The three Unes show paraboUc fits over intervals of 150, 40 and 
5 frames. 



sweeping time Tg) we have observed that the slope of the ao vs. r/ curve has some dependence 
on the finite size of the measurement volume, although the extrapolation to r/ = is found to 
be independent of the measurement volume. For this reason the r/ — > extrapolation is a more 
reliable measure of ao than its value at any finite r/. 

The variance of the (unnormalized) x component of the acceleration is shown as a function of 
the rms velocity u in figure 25. It is found that the predicted scaling 

= (4.8) 

(From equation 4.4, using e = /L), is observed over nearly 7 orders of magnitude in accel- 
eration variance, or nearly 2 orders of magnitude in velocity variance. The scahng of ao with 
Rx is plotted in figure 26 and (tabulated in Table 4.2). In this plot, a constant value of ao would 
indicate Heisenberg-Yaglom scaling. As Rx is decreased below 500 the value of ao decreases 
substantially showing a small but significant departure from universal scaling. The dependence 
of ao on Rx in this regime seems to be quahtatively consistent with DNS results. It is not sur- 
prising that DNS and experimental results do not match exactly at small Reynolds number, since 
in this range the acceleration is coupled to the large scales of the flow. The DNS results were ob- 
tained for isotropic turbulence with periodic boundary conditions, which differs markedly from 
the anisotropically forced turbulence between counter-rotating disks used in the experiment. At 
high Reynolds number {Rx > 500) ao appears to be independent of Rx, which is consistent with 
Heisenberg-Yaglom scahng. Due to experimental uncertainties, very weak deviations such as the 
Rx'^'^^^ prediction of the Borgas multi-fractal model cannot be ruled out by this data. 

Please note that it has been shown that the mean squared pressure gradient in a turbulent 
flow — and therefore the mean squared acceleration, and ao — is closely related to the fourth 
order velocity structure functions and to the inertial range flatness factor (Hill & Wilczak 1995). 
This would indicate that the ao vs Rx curve in figure 26 could be compared with the flatness 
factor (F) vs Rx curve reported by Belin et. al. in another experiment involving a flow between 
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Figure 23. Normalized acceleration variance (ao) of x component as a function of fit interval normalized 
by the Kolmogorov time (t//t,,). The circles represent data taken at R\ = 970, the solid line indicates the 
best fit of the function in equation 4.7, and the dashed line shows the exponential term. A lOextrapolation 
to zero fit time: 
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Figure 24. Plots of the normalized variance of an acceleration component, ao, as a function of nor- 
malized fit time T//r^ for the full range of Reynolds number. The lines indicate the exponential term 
C exp Dt + Et^ from the best fit for each curve. The tabulated values of oo are derived from the extrapo- 
lation of this term to r = 0. (a) x component (b) z component. 
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Figure 25. The variance of the acceleration plotted as a function of w. The straight line is the best fit of 

M^/^ (g3/2 with e = M^/L) to the data. 

counter-rotating disks (Belin et at. 1997). Some similarity is evident, despite the fact that the 
flow configurations are quite different (Hill 2000). 

It is also evident that larger qq values are obtained for the transverse component than for the 
axial component of the acceleration, and that the level of anisotropy decreases as the Reynolds 
number is increased. The anisotropy of the acceleration variance is best illustrated by calculating 
the ratio {ao)x/(ao)z as a function of Rx, as shown in figure 27. At the lowest Reynolds number 
of 140, (ao)x/ {(io)v = 1-52, indicating that the ratio of the standard deviations of the acceler- 
ation components is 1.23. As the Reynolds number is increased to 970, the level of anisotropy 
decreases to a very small value. This observation is consistent with recent experimental results 
which indicate that anisotropy persists to high Reynolds numbers (Kurien & Sreenivasan 2000; 
Shen&Warhaft2000) 

It is also informative to compare figure 27 with the inset to figure 14, which shows the ratio of 
the standard deviations of the velocity components. On the one hand, the particle velocities are 
associated mainly with the large scales and the level of anisotropy is large and independent of 
Reynolds number, as expected. On the other hand, the particle accelerations come mainly from 
the dissipation range scales. As expected, the level of anisotropy for the acceleration is smaller 
than for the velocity, and decreases as the Reynolds number is increased. 

4.3. Acceleration Flatness 

It is also interesting to quantify the degree of intermittency of the acceleration by calculating the 
acceleration flatness ((«|) / (af )^) as a function of Reynolds number. In this case, one encounters 
the same difficulty as for the variance; the flatness varies as a function of the interval over which 
the acceleration is calculated. One must also confront the added difficulty that the tails of the 
acceleration distribution are so long that prohibitively large data samples would be needed to 
definitively converge the fourth moment, as illustrated in figure 20 above. In view of figure 20, it 
seems possible that the unconverged tails of the acceleration distribution could make substantial 
contributions to the fourth moment. In addition, the ability of the tracer particles to fully follow 
rare violent events which contribute most to the flatness has not been fully estabhshed. Therefore, 
we can only set lower bounds to the flatness. 

The flatness for the transverse component of the acceleration is shown as a function of fit time 
in figure 28. An increasing trend is evident as r/ decreases, but the dependence on t/ is more 
difficult to define than in the analogous curves for the acceleration variance (figure 24). We have 
not found it possible to make a formal extrapolation for t/ ^ in this case, but by tabulating 
the flatness at a value of r/ shghtly above the onset of power-law position uncertainty we can 
compile a plot of lower bounds on the acceleration flatness, shown in figure 29. The flatness is 
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Figure 26. The Kolmogorov constant oo calculated for transverse (filled circles) and axial (filled squares) 
components of the acceleration, as a function of Reynolds number. These values have a 10% correction 
applied to them as discussed in Section 4.5. Values obtained from direct numerical simulation of turbulence 
by Vedula and Yeung(Vedula & Yeung 1999), and by Gotoh and Rogallo(Gotoh & Fulcayama 2001) are 
also shown by open squares and open triangles, respectively. 




found to be at least 25 at the lowest Reynolds numbers studied and the lower bound increases to 
approximately 60 as the Reynolds number is increased. 

These large flatness values indicate that the acceleration is more intermittent than the other 
small scale quantities in turbulence. For example, at R\ = 200, the longitudinal velocity deriva- 
tive flatness is 6.0 (Vanatta & Antonia 1980), and the scalar gradient flatness is 17 (Shen & 
Warhaft 2000), while the acceleration flatness in both experiments and simulation (Vedula & 
Yeung 1999) is approximately 30. 

4.4. Acceleration Correlations 

Having measured the acceleration as a function of time for the particle tracks, it is straightforward 
to evaluate the acceleration autocorrelation function, 

where af{t) is the acceleration component along a particle trajectory and () denotes averaging 
over t for an ensemble of tracks. In principle, the denominator can be simplified and expressed 
as the variance of af but because the lengths of the tracks are of the same order as the range of 
T, end effects can be significant; (a^(t)^) and {af{t + r)^) are compiled from different data 
samples and can differ sUghtly, particularly at large r. 

The acceleration autocorrelation is shown for R\ = 970 and R\ = 690 in figure 30. In the 
main graph the time axis is normalized by the Kolmogorov time t,, for each curve, and in the inset 
the unnormalized time axis is in units of seconds. The Kolmogorov time is proportional to u~^^^ 
and changes by a factor of 2\/2 between the two values of R\ studied. The data reproduces 
this scaling to an accuracy of a few percent. Calculations of the acceleration autocorrelation 
function in direct numerical simulation of turbulence at Rx = 140 indicate that it crosses zero at 
2.2r,, (Yeung 1997). In our experiment we see the zero crossing at 2.1r^ at R\ = 970 and 3r^ at 
R\ = 690. The data collapse shown in figure 30 may be regarded as an independent confirmation 
of the measurement of the energy dissipation, described above in Section 3.3. 

The autocorrelation data is compiled using the same trajectories as were used for the study of 
the acceleration moments. This presents some difficulty because of the small size of the mea- 
surement volume. The autocorrelation function is compiled under the impficit assumption that 
the particle remains in view after a time interval r has elapsed, however, particles with large ac- 
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f Rx 
(Hz) 



(ao)x 



{ao)z 



{ao)x/{ao)z 



0.15 140 4.60 ±0.40 3.01 ± 0.27 

0.30 200 4.99 ±0.28 4.07 ± 0.48 

0.60 285 5.26 ±0.28 4.49 ± 0.70 

1.75 485 6.01 ±0.33 5.12 ±0.23 

3.50 690 6.18 ±0.42 5.85 ± 0.27 

7.00 970 5.61 ±0.35 5.23 ± 0.34 



1.53 
1.23 
1.17 
1.17 
1.06 
1.07 



Table 3. ao as a function of Reynolds number for transverse and axial acceleration. 



celerations are less likely to remain in view than particles with small accelerations. The ensemble 
of measurements contributing to autocorrelation function at r = 5t,, have a standard deviation 
that is a factor of 0.7 smaller than that of the unconditional acceleration variance, indicating that 
the bias towards low acceleration events is appreciable. To make a definitive study of the acceler- 
ation autocorrelation, and particularly to obtain results at low Reynolds number, a more powerful 
laser which can illuminate a larger measurement volume will be required. 



As discussed before, measurements of Lagrangian statistics are comphcated by the fact that all 
statistics are implicitly conditional on the event being successfully measured. For measurements 
of acceleration, the particle must remain in view for a sufficient time for its acceleration to be 
determined, which may introduce biases due to measurement volume effects. The extrapolation 
to zero fit time further complicates the matter, because there are two effects which may affect 
the measured acceleration variance as the fit interval is varied. As fit interval is decreased, the 
accelerations measured for particular trajectories will increase because there is less coarse grain- 
ing over turbulent fluctuations and noise. In addition, as the fit interval is decreased the sampling 
of trajectories becomes inclusive because the particles need not remain in view for as long an 
interval. 

To determine how the extrapolation to zero fit time is related to the true acceleration variance, 
we used a simulation of the detection process. The simulation is based on three dimensional 
particle trajectory data from Vedula and Yeung's DNS simulation (Vedula & Yeung 1999) at 
Rx = 240, which is known to have ao = 3.44. These trajectories were dimensionalized using the 
viscosity of water and the rate of energy dissipation which would be needed to produce the same 
Reynolds number in our apparatus. These scaled trajectories were used as input to a computer 
model which simulated the illumination of the particles and the imaging of our measurement 
volume onto the strip detectors, including diffraction and defocusing of the imaging system. 
The computer model also simulated charge collection by the strip detector, including correlated 
and uncorrelated noise, charge diffusion, and inoperative pixels. The intensities produced by the 
simulation were then run through the real-time thresholding and compression algorithms used in 
the experiment, and were subsequently processed by the same data processing algorithms used 
for actual experimental data. 

The measurement of ao as a function of normalized fit interval is shown in figure 3 1 (which 
may be compared with figure 23, above). The exponential dependence on r/ is reproduced for 
large r/, as is the power law dependence for small r/. The extrapolation for r/ — > gives 
a value of 3.80. This value is 10% higher than the value of 3.44 which was measured directly 
from the simulation (Vedula & Yeung 1999). All normalized acceleration variance measurements 
presented in this paper have been rescaled to correct for this overestimation. 



4.5. Simulation of Detection process. 
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Figure 27. The ratio of ao calculated for transverse {x) and axial {z) components of the acceleration as a 

fimction of Reynolds number. 

5. Effect of Finite Particle Size on Particle Accelerations 

An important question which must be considered is the extent to which the polystyrene tracer 
particles are equivalent to ideal fluid particles. In order to address this question, we have repeated 
our acceleration measurements using tracer particles with a range of particle diameters and fluid 
densities. Not only are these measurements essential to validate the Lagrangian measurements, 
but the data for large particle sizes offers a new perspective on the motion of finite size particles 
in turbulence. 

Phenomenologically, we can identify two mechanisms by which a discrepancy between statis- 
tics of tracer particle trajectories and fluid particle trajectories could arise. (1) It is possible that 
the tracer particle's density is different from that of the fluid, so that it experiences an acceler- 
ation that is different than the acceleration that would be experienced by a fluid particle in the 
same location. We call this the density dependent effect. (2) It is possible that the tracer particle's 
size is large enough that the unperturbed flow would change significantly over the volume of the 
particle. Since this effect remains even when the tracers are density matched, we call this the 
density independent effect. The dimensionless numbers that measure these effects are the ratio 
of the particle to fluid density, Pp/p/, and the ratio of the particle diameter to the smallest length 
scale in the flow, d/ij. A commonly used parameter is the ratio of the Stokes time {(Ppp/vp /) to 
the Kolmogorov time {rj^ /v) which is the combination of the size and density parameters. 

Particle motion in fluid flows has been extensively studied both as a fundamental fluid dynam- 
ics topic (Basset 1888; Corrsin & Lumley 1956; Maxey & Riley 1983) and in order to vaUdate 
measurement techniques that rely on tracer particles (Buchave et al. 1979). Most of the work in 
this area has focused on the dynamics of particles when the density independent effect is neghgi- 
ble. Not only is this regime amenable to exact theoretical modelling, but it also is a good approx- 
imation for the tracers used in many particle image velocimetry and laser Doppler anemometry 
experiments. In these experiments, it is desirable to have the spacing between particles less than 
or equal to the size of the smallest flow structures, so the particle size must be much smaller than 
the smallest structures. 

In order to determine the effects of finite particle size on the acceleration measurements, we 
repeated the measurements at R\ = 970 using particles of diameters between 26 /xm and 450 /itm 
( 1.44 T] and 26 rf). The Stokes times for these particles ranges from 0.67 ms to 216 ms, which 
should be compared to the Kolmogorov time of 0.32 ms. We also measured acceleration of 
450 jum particles in NaCl solutions with densities between 1 .00 g cm~^ and 1 . 1 1 g cm~^ . Since 
the particle density is 1 .06 g cm~^, this range of densities encompasses both negative and positive 
particle buoyancy. 

In figure 32, the measured Kolmogorov constant ao is shown as a function of particle size for 
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Figure 28. Flatness of the transverse component of the acceleration as a function of fit interval t/. 
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Figure 29. Lower bound on the flatness of the transverse and axial components of the acceleration as a 
function of Rx. The error bars indicates an estimated uncertainty due to the inability to extrapolate t/ to 
zero, but do not reflect uncertainty in the convergence of the fourth moment. 
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Figure 30. The acceleration autocorrelation function is shown for data at Rx of 970 and 690, where 
the time axis has been scaled with the Kolmogorov time r,,. The accelerations are calculated by direct 
differencing of trajectories which have been smoothed with a Gaussian kernel, where the standard deviation 
of the kernel is 0.5 r,,. The insert shows the same data in which the time is unsealed and plotted in seconds. 
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Figure 3 1 . The plot shows ao as a function of normahzed fit time r/ /r,, for trajectories obtained from 
direct numerical simulation of turbulence(Vedula & Yeung 1999). The value of ao obtained from direct cal- 
culation from the simulation is 3.44. The extrapolation technique overestimates this value by approximately 



acceleration data taken with pure water at R\ — 970. A 70% decrease is seen in the measured 
acceleration variance between the smallest and largest particles used. 

The dependence of the particle acceleration on the density mismatch is shown in figure 33, 
where oq is plotted as a function of relative density for 450 /im diameter particles at R\ = 970. 
In this case, the density of the fluid is varied by the addition of NaCl, and the values of ao 
are corrected for the slight change in kinematic viscosity which occurs as the density is varied. 
The data for pure water appears at effective particle density pcS = 1.06, and reflects the strong 
underestimate of particle accelerations with large dense particles. If the suppression of ao were 
due to density mismatch, the effect would reverse in the case where the particles are lighter than 
the fluid they displace, and the ao values would exceed the small particle limit. The data at poff = 
0.96 in figure 33 (taken at high salt concentration) shows that although aq increases somewhat, it 
still remains far below the small-particle limit of 6 as seen in figure 32. Figure 34 shows the PDF 
of the acceleration for different particle sizes and fluid densities. Out to 10 standard deviations, 
the shape of the distribution is not strongly affected by particle size. This data demonstrates that, 
in our experiment, the effective coarse graining of the acceleration by the finite size particles is 
more important than the differential buoyancy due to the density mismatch in our measurements. 

A simple model based on the K41 phenomenology can be used to predict the scaling of the 
acceleration with particle size in this regime where density mismatch is not important. We as- 
sume that the effect of finite particle size is that only flow structures larger than the particle size 
contribute to the acceleration. This means that the scaling of the acceleration variance can be 
determined by starting with, < a^ >~ e'^/^j/^^/^, and replacing the viscosity with the value that 
would make the Kolmogorov scale equal to the particle diameter, ly e^/'^d^/'^. Consequently, 
the acceleration variance should scale as d^^^'^ for large particle sizes. We expect that below 
some particle size this scaling is no longer adhered to, and the acceleration becomes independent 
of particle size. The solid line in figure 32 shows the d^^^^ scaling. There is not enough data to 
confirm this theory, but the data is consistent with this scaling. For large particle size, the data 
can be interpreted as following the d^^^^ scaling, and for particle size less than 100 /im (5 77) 
there is a turn over to the acceleration being independent of particle size. 

To provide an upper limit on the possible deviation of tracer particle acceleration from that 
of fluid particles, we note that the dependence of ao on particle diameter in figure 32 could be 
interpreted as being roughly linear. An extrapolation of the linear dependence to zero particle 
size indicates that the acceleration variance of the 46 /um particles used for the acceleration 
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Figure 32. Normalized acceleration variance (oo) as a function of particle diameter at R\ = 970. Open 
circles with error bars are for particles of relative density 1.06. The closed circle with error bars is for 
particles of relative density 1.0. The Kolmogorov length is 18 /^m and the relative density is defined as 
peS = pparticie/pfluid- The solid line shows scaling. 

measurements reported above is within 7% of that of fluid particles. This is really a worst case 
scenario, since the values obtained for 26 /xm and 46 /xm particles are indistinguishable within 
experimental uncertainties, and it is expected that the particle size dependence has zero slope at 
the origin. Note that this data is at the largest Reynolds number reported in this paper, and the 
particle size dependence will be weaker at lower Reynolds numbers. We conclude that the 46 /xm 
particles are small enough to measure the acceleration variance at the 10% accuracy which we 
have reported. 

6. Conclusion 

We have reported high resolution measurements of the Lagrangian properties of a high Reynolds 
number turbulent flow. The particle tracking system, based on the CLEO III silicon strip detector, 
is able to measure positions with spatial resolution of 0.5 yum (1/40 of the Kolmogorov distance, 
or 1/6400 of the field of view) and temporal resolution of 14 fxs (1/20 of the Kolmogorov time) 
in an Rx = 970 flow. 

The primary results of this paper concern the fluid particle accelerations which, unlike the 
particle velocity, cannot be measured with an Eulerian (fixed) probe. We tested the long-standing 
Heisenberg-Yaglom prediction that the acceleration variance scales as e^/^, (or u^/^). At low 
Reynolds number, our data is consistent with DNS. However, we find that for 500 > Rx > 970 
the value of ao is approximately constant, indicating that K41 theory approximates the acceler- 
ation scaling in this range of Rx well. Turbulent intermittency models generally predict a very 
weak dependence of ao on Reynolds number (for example, R^-^"^^ in Borgas' model) but exper- 
imental errors currently make it impossible to test such predictions. 

The fact that the flow between counter-rotating disks has a well-defined anisotropy allows us 
to test the coupling of the acceleration to the structure of the large scale flow. While the veloc- 
ity exhibits a significant anisotropy which is independent of Reynolds number, the acceleration 
exhibits a smaller anisotropy which decreases as the Reynolds number is increased, as shown in 
figure 27. 

The autocorrelation function of the acceleration component was also measured. Because of 
measurement volume effects, this measurement was restricted to relatively high Reynolds num- 
bers. The functions obtained agree with low Rx DNS simulations, crossing zero at approximately 
2.1 T^, and exhibit Kolmogorov scaUng of the time. 

The acceleration is found to be an extremely intermittent variable, as evidenced by the very 
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Figure 33. oo (normalized acceleration variance) as a function of relative particle density for 450 /um 

particles at Rx = 970. 



long stretched exponential tails of the acceleration component PDF, shown in figures 19 and 21. 
The level of intermittency, reflected in the acceleration component flatness, is found to increase 
with Reynolds number. This is consistent both with the general expectation that turbulent inter- 
mittency should increase as the width of the scaling range increases, and with DNS measurements 
of the flatness at low Reynolds number. 

The fast temporal response of the particle tracking system allows us to resolve the extremely 
violent events that make up the tails of the acceleration component PDF. The trajectory in figure 1 
represents one such event, in which the particle acceleration rises to 16,000 ms~^, 1,600 times 
the acceleration of gravity, and 40 times the rms value. Although we do not have quantitative 
evidence, we make the observation that the high acceleration events appearing in the tails of 
the acceleration distributions (such as figure 1) seem to be associated with coherent structures 
which persist for many Kolmogorov times, substantially longer than the correlation time of the 
acceleration components. 

Particular attention has been paid to the determination of measurement errors. We have devel- 
oped a numerical simulation of the detection process that takes input trajectories and models the 
optics and electronics to create artificial data. Using this approach we have measured the extrap- 
olation error in the acceleration variance data and corrected for it. We have also confirmed that 
our analysis codes and choices for the analysis parameters are not biasing the results. 

The dependence of the particle acceleration variance on the tracer particle size and density 
was also measured. It is found that the acceleration observed with 46 jim diameter polystyrene 
particles is within 7% of the value that would be obtained from ideal fluid particles, even at the 
highest Reynolds numbers studied. We find that particles of diameter 450 jim (26 rf) have an 
acceleration variance that is a factor of 3.6 smaller than fluid particles. In our experiment, this is 
caused primarily by the size of the particles and is only sUghtly affected by their density mismatch 
with the fluid. These measurements validate our Lagrangian acceleration measurements, and 
in addition they highhght the need for a deeper understanding of the motion of large density 
matched particles in turbulence. 

The techniques and measurements presented in this paper suggest many possibilities for the 
future development of optical particle tracking in turbulence. The use of high resolution imaging 
equipment and careful attention to measurement errors will continue to be essential in future mea- 
surements. Possible extensions of the techniques presented here include lengthening the tracking 
times, improving the resolution of the three-dimensional measurements and tracking larger num- 
bers of particles simultaneously. These will allow precise high Reynolds number experimental 



32 Greg A. Voth, A. La Porta, Alice M. Crawford, Jim Alexander, Eberhard Bodenschatz 











• 1.105 ms 1 


a « 


o 0.005 ms 




A 0.006 ms 


/\ 




1 \ 

























a/(a2)l/2 



Figure 34. PDF of the acceleration for different particle sizes and fluid densities at R\ = 970 normalized 
by the rms acceleration. Closed symbols are for particles of size 50 /irn and open symbols are for particles 
of size 450 /irri. Circles are for peff = 1.06. Triangles are for poff = 0.96. The unnormalized standard 
deviation for each curve is shown in the graph. Out to 10 standard deviations the shape of the curve is not 
strongly affected. 



measurements of additional quantities including relative dispersion, scaling in the inertial time 
range, and the geometry of multi-particle Lagrangian motion. 
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Appendix A. Extraction of Particle Tracks from Intensity Data 

The raw data from the strip detector consists of a series of one-dimensional intensity profiles 
taken at regular time intervals. The task of recognizing peaks in these intensity profiles and as- 
sembling them into two-dimensional or three-dimensional tracks is, in principal, similar to the 
case of two-dimensional CCD images (Virant & Dracos 1997; Ott & Mann 2000; Dracos 1996). 
However, the strip detector poses unique challenges. Particle tracks cross much more often in 
one-dimensional than in two-dimensional, and flaws in the detector lead to a number of inop- 
erative strips. Both of these effects cause frequent drop-outs in tracks, which would result in 
fragmentation of the trajectories. 

The strategy for extracting tracks from the intensity data is as follows. The task is divided 
into four phases, peak detection, track building, track splicing and track filtering. In cases where 
two-dimensional trajectories are investigated, the tracks are then passed to a matching algorithm 
which associates x and y tracks to form two-dimensional trajectories. 

Peak Detection. The algorithm scans each frame and identifies distinct groups of above-threshold 

strips, and searches each group for one or more peaks, taking into account the existence of in- 
operative strips. Vahd peaks are passed to a routine which calculates moments of the peak and 
fits it to a Gaussian function, returning properties including peak centre, amplitude, full-width at 
half maximum, area skewness and flatness. The output of the peak detection algorithm is a Ust of 
peaks for each frame. 

Track Building. The track finding algorithm assembles the peaks found by the peak detector 
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into time-continuous tracks. This algorithm builds tracks incrementally, extrapolating each track 
forward one frame and searching for a peak which is the continuation of the track. If any ambi- 
guity is found — if there is more than one peak which might continue the track, or if there is more 
than one track which might continue onto a peak — ^the track building algorithm fails to resolve 
the match and begins new tracks as necessary. 

Splicing. The list of tracks generated by the track builder is passed to a splicer which connects 
track fragments which the track builder failed to connect due to ambiguities. Such ambiguities are 
quite common in one-dimensional projection, since tracks can cross in one-dimensional projection 
even when the two-dimensional trajectories are distinct. The track builder looks more than one 
frame ahead, and compares forward extrapolations of the ends of tracks to backward extrapola- 
tions of the beginnings of tracks. It uses an iterative algorithm to make the best matches, while 
leaving ambiguous matches unresolved. 

Filtering. After the spUcing algorithm has exited, the tracks are processed by a filtering algo- 
rithm. The algorithm will delete data points which meet certain criterion. These criterion include 
upper and lower limits on the amplitude, and upper limits on the width and flatness of the peak. 
The algorithm also deletes data points which coincide with a list of inoperative strips. The idea 
is to keep marginally accurate data during the track assembly process to maintain the continu- 
ity of the tracks, but to exclude such data from subsequent statistical analysis, where the large 
uncertainties could be detrimental to the analysis. 

Coordinate Matching. Most of the acceleration data reported below is derived from two- 
dimensional trajectories. To obtain two-dimensional trajectories, it is necessary to match the x 
vs. t and z vs. t trajectories which are obtained from the two strip detectors. All such data is 
taken in the configuration shown in figure 5(a), in which the same image is projected on two strip 
detectors, and only the segments of tracks thai lie in the overlap region of the two detectors (See 
figure 7) are used. In this case, the optical intensity recorded by the two detectors for a single 
track is highly correlated. The autocorrelation function between the two trajectories is therefore 
a good figure of merit for matching the tracks. 

Appendix B, Reinterpretation of previous data. 

The data presented in this paper should be compared with previously pubhshed measurements 
of the same quantities in the same turbulent flow, but using a conventional position sensitive 
photodiode instead of the strip detector to measure particle trajectories (Voth et al. 1998). The 
previous study reached the correct conclusion that the Heisenberg-Yaglom scaling of the accel- 
eration variance is observed at high Reynolds numbers, but because of limitations of the detector 
technology the numerical values obtained for the physical quantities involved were inaccurate. 

The discrepancy can be largely attributed to the high noise level in the position sensitive photo- 
diode that resulted in much greater uncertainty in the particle position measurements and which 
necessitated the use of much larger tracer particles. The intrinsic bandwidth of the UDT DLS-10 
position sensitive photodiode used in the previous experiments is approximately 5 MHz, but tem- 
poral averaging of the signal was used to increase the effective position resolution at the expense 
of temporal resolution. Signal averaging with an effective bandwidth of 100 kHz was used in 
conjunction with tracer particles of diameter 250 ^m or 450 ^m, which gave a position resolu- 
tion of 10 /Ltm over a field of 1.5 mm, corresponding to a dynamic range of 150. By contrast, the 
strip detector at 70 kHz has a position uncertainty of approximately 0.1 strips on a field of 512 
strips (0.8 /Ltm over a field of 4 mm), giving a dynamic range of 5000, which is a factor of 30 
improvement. Another critical advantage of the strip detector is that it is able make these posi- 
tion measurements with 46 /im diameter particles. As reported in Section 5 above, the 250 /im 
particles which were required by the position sensitive photodiode do not follow the flow at high 
Reynolds number and result in significant underestimate of the acceleration. 
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The position sensitive photodiode is strictly limited to single particle statistics, so that the en- 
ergy dissipation could not be measured in terms of velocity structure functions. Instead, the ac- 
celeration autocorrelation function was compared with DNS results to estimate the Kolmogorov 
time, and this value was used to determine the rate of energy dissipation. Presumably due to 
the fact that the 250 jim particles did not adequately follow the flow, the estimate of the Kol- 
mogorov time was a factor 2.37 larger than the current value, and using e = (j//e)^/^, this gave 
an estimate of the energy dissipation which is a factor 5.63 below the current value at a given 
propeller speed. The rms velocity u was also underestimated by 10% in the old experiment, and 
using equation 3.7 the new values obtained for Rx are smaller than the old by a factor of 0.52, 
so the previously reported range of Rx (985-2021) would correspond to (512-1037) using the 
new calibration. The ability to work at much smaller Rx (down to 140) results from the use of 
smaller particles, which do not settle out of the flow even at very slow propeUer speeds. 

The Kolmogorov constant oq is calculated from 




(Bl) 



and is also affected by this inaccuracy in the estimate of e, and would have increased the value 
of ao by a factor of 17. However, the particle size effects and low spatial resolution caused the 
acceleration variance to be underestimated by a large factor. In the old study, an extrapolation 
to zero fit time was attempted, but large position errors obscured the exponential dependence of 
acceleration variance on fit time that is seen in figure 23. What was interpreted as "measurement 
noise" was in reality a blend of measurement error and and short-time turbulence contributions. 
The algorithm used to calculate the acceleration variance gave a value which was effectively 
coarse grained over a time interval which was beUeved to be 2.5 using the old calibration, but 
which is actually 6 using the new calibration. As may be seen from figure 23, this results in 
a substantial underestimate of the acceleration variance. The acceleration variance was further 
reduced by the particle size effect associated with the use of 250 /zm particles. The fact that the 
value of ao reported in the previous experiment (7 ± 3) is close to the value now reported (« 6) 
is due to the fortuitous cancellation of errors in the measurement of the dissipation rate and of 
the acceleration variance. 
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